UNITED  STATES  AIR  FORCE 
AIR  UNIVERSITY 

AIR  FORCE  INSTITUTE  OF  TECHNOLOGY 

Wright«Patt*rson  Air  Fore*  Bat*, Ohio 


DISTRIBUTION  STATEMENT  A 
Approved  tot  public  release; 
Distribution  Unlimited 


80  3  14  027 


AF IT/GAE/AA/7  8D- 1 1 


0  LEVEL 


SLIP  DAMPING 
OF 

TURBINE  BLADES 

THESIS 

AF IT/ GAE/AA/7 8D- 1 1 v  Michael  C.  Kimberling 

Capt  USAF 


DISTRIBUTION  STATEMENT  A 

Approved  lor  public  release; 
Distribution  Unlimited 


DTIC 


1 


Acknowledgements 

I  would  like  to  take  this  opportunity  to  express  my  gratitude  to 
some  of  those  individuals  who  were  most  prominent  in  assisting  me  to 
complete  this  effort.  My  thesis  advisor,  Dr.  Peter  J.  Torvik,  who  as¬ 
sisted  me  in  understanding  many  aspects  of  this  study,  and  Dr.  Clarence 
Edstrom,  who  assisted  me  with  much  of  the  mathematical  work  done  for 
this  study,  are  both  deserving  of  recognition  for  their  efforts.  The 
same  must  be  said  for  Captain  Dan  Hayes  who  was  invaluable  in  helping 
me  with  the  experimental  portion  of  this  study.  Also  worthy  of  note 
was  the  excellent  support  provided  by  the  AFIT  shops  in  the  manufac¬ 
turing  of  the  experimental  model. 


ACCESSION  tor 

NTIS 

DOC 


Whit*  Section 
Buff  Suction 


UNANNOUNCED 

JUSTIFICATION 


^ISII«UTW^W«lin  CUB 


r  i  i  1 


□  □ 


Contents 


Page 


Acknowledgements  .  ii 

List  of  Figures . v 

List  of  Tables . vi 

List  of  Symbols . vii 

Abstract .  ix 

I.  Introduction  .  1 

II.  Experimental  Development  .  3 

Model  Description  .  3 

Test  Set-up .  8 

Test  Procedure .  8 

Data  and  Discussion .  11 

III.  Analytical  Development  .  18 

Model  Description .  18 

Contact  Surface  State  of  Stress  .  18 

Contact  Surface  Displacements  .  23 

Points  of  Slip  Initiation,  and  Zero  Normal  Force  .  .  31 

Energy  Dissipated  by  Slip .  33 

Energy  Dissipated  by  Air  Damping  .  33 

Energy  Dissipated  by  Hysteresis  .  34 

Restoring  Moment  Due  to  Slip .  35 

Equation  of  Motion  and  Boundary  Conditions  .  36 

Solution  to  Equation  of  Motion  .  37 

IV.  Recommendations .  44 

Bibliography  .  45 


iii 


Figure 


Page 


List  of  Figures 

1  Picture  of  Experimental  Blade  Model  .  4 

2  Dimensions  and  Sign  Convention  of  Both  Models  .  5 

3  Picture  of  Experimental  Disk  Model  .  4 

4  Picture  of  Assembled  Disk-Blade  Model  (side  view)  ....  6 

5  Picture  of  Assembled  Disk-Blade  Model  (end  view)  ....  6 

6  Picture  of  Spring  Beam .  7 

7  Schematic  of  Test  Set-up  .  9 

8  Picture  of  Test  Set-up  in  Lab .  10 

9  Plot  of  Force  Versus  Frequency  for  P  =  10  lbs  and 

Variations  in  Tip  Displacement  .  12 

10  Plot  of  Force  Versus  Frequency  for  P  =20  lbs  and 

Variations  in  Tip  Displacement  .  14 

11  Plot  of  Force  Versus  Frequency  for  Displacement  = 

.0001  in.  and  Variations  in  P .  16 

12  Plot  of  Force  Versus  Frequency  for  Displacement  = 

.0003  in.  and  Variations  in  P .  17 

13  Breakdown  of  Applied  Loads  .  19 

14  Typical  Stress  Distribution  of  Contact  Surface  .  24 

15  Relative  Displacement  Across  Contact  Region  .  30 

16  Schematic  of  Calibration  Set-up  for  the 

Accelerometer  and  Impedance  Head  Calibration  .  50 

17  Picture  of  Force  Table  and  Guide .  52 

18  Picture  of  Force  Table  in  Place  on  the  Disk  Model  ....  52 


v 


List  of  Tables 


Table  Page 

1  Listing  of  Eigenvalues  and  Natural  Frequencies 

For  the  Homogeneous  Beam .  48 


vi 


List  of  Symbols 


P 

0 

L 

F 

r 

t 

b 


1>  v2»h>3 

Q,  R 


r0 

M 


<ST 

<5U 

<$WK 

Y 

I 

o 

W(X,T) 

E 

I 


w,ii2 


S) 

v 


Force  applied  by  setscrew  and  spring  beam  (lbf) 

Angle  defining  the  contact  surface  of  the  disk  model 
Length  of  blade  (in) 

Magnitude  of  sinusoidal  blade  tip  force  (lbf) 

Radius  of  blaae  root  (in) 

Depth  of  blade  root  (in) 

Width  of  blade  model  (in) 

Stress  functions 

Points  of  slip  initiation  on  contact  surface  (deg) 

2 

Radial  stress  (lb/in  ) 

2 

Tangential  stress  (lb/in  ) 

2 

Shear  stress  (lb/in  ) 

Restoring  moment  due  to  slip  (in-lb) 

Dynamic  coefficient  of  friction 
Static  coefficient  of  friction 
Variation  in  kinetic  energy  of  blade 
Variation  in  potential  energy  of  blade 
Variation  in  work 

Mass  per  unit  length  of  blade  (lbm/in) 

2 

Mass  moment  of  inertia  of  blade  root  about  point  0  (lbm-in  ) 
Displacement  of  the  blade  in  the  y  direction 

2 

Modulus  of  elasticity  of  the  blade  (3861bm/sec  in) 

4 

Area  moment  of  inertia  of  the  blade  (in  ) 

Forcing  frequency 

Natural  frequency  of  system 


Poisson's  Ratio 


AFIT/GAE/AA/78D- 11 


Abstract 

An  experimental,  and  an  analytical,  model  of  Turbine/Compressor  blade 
with  a  circular  root  geometry  is  developed  for  the  purpose  of  studying 
the  effects  of  slip  damping  on  the  response  of  the  blade.  The  expression 
for  the  energy  dissipated  by  slip  damping  is  presented.  In  addition  the 
expressions  for  the  energy  dissipated  by  air  and  hysteresis  damping  are 
presented.  The  equation  of  motion  for  the  blade  is  developed,  and  solved, 
to  provide  an  expression  relating  the  various  parameters  of  the  experi¬ 
mental  and  analytical  models. 


IX 


I.  Introduction 


Today  there  are  significant  problems  associated  with  the  failure 
of  Turbine/Compressor  blades.  The  failures  are  due  to  increasing  de¬ 
mands  in  the  areas  of  performance  and  service  life  from  the  jet  engines. 
To  increase  the  performance,  higher  thrust  to  weight  ratios  are  deman¬ 
ded.  As  a  result,  the  excess  material  of  the  Turbine/Compressor  blades 
is  being  reduced  to  the  absolute  minimum  and  the  margins  of  safety  are 
being  reduced  to  the  minimum  allowable  limits.  Longer  service  life  re¬ 
quirements  have  brought  about  problems  with  fatigue  which  are  difficult 
to  predict  due  to  the  random  nature  of  the  dynamic  loading. 

The  environment  of  the  Turbine/Compressor  blade  is  very  severe. 

The  stresses  generated  within  the  Turbine/Compressor  blades  are  very 
high.  The  loading  which  acts  on  the  blades  is  very  dynamic  and  random. 
The  loading  is  expecially  severe  at  the  resonant  frequencies  of  the 
Turbine/Compressor  blades  where  there  is  increased  displacement  which 
results  in  an  increase  in  the  stresses.  In  addition,  for  Turbine  blades, 
there  is  the  additional  problem  of  a  very  high  temperature  environment. 

It  is  important  to  understand  the  mechanisms  which  operate  to  reduce 
the  effect  of  this  dynamic  loading.  With  a  knowledge  of  these  damping 
mechanisms,  it  might  be  possible  to  exploit  them  to  obtain  a  further 
reduction  in  the  effect  of  the  dynamic  loading. 

There  are  three  damping  mechanisms  which  operate  to  reduce  the  re¬ 
sponse  of  a  Turbine/Compressor  blade.  They  are  hysteresis  or  material 
damping,  air  damping,  and  slip  damping.  This  paper  is  concerned  with 
slip  damping  and  so  the  other  two  forms  of  damping  were  addressed  only 
to  the  degree  necessary  to  provide  the  means  of  isolating  the  effect  of 
the  slip  damping. 
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This  study  was  undertaken  to  establish  both  an  experimental,  and  an 
analytical  model  of  a  Turbine/Compressor  blade.  These  models  may  be  used 
to  study  the  effect  of  slip  damping  on  the  response  of  a  system  similar 
to  a  Turbine/Compressor  blade  mounted  in  a  disk.  The  analytical  and 
experimental  models  might  also  be  used  to  study  ways  of  increasing  the 
slip  damping  which  occurs  on  the  contact  surface  between  the  blade  root 
and  the  disk.  For  example,  if  the  slip  damping  mechanism  could  be  main¬ 
tained  at  higher  engine  RPM  levels  there  would  be  a  significant  contri¬ 
bution  to  the  reduction  of  the  response  of  the  Turbine/Compressor  blade. 

The  primary  area  of  investigation  presented  in  this  paper  is  the 
development  of  the  analytical  model.  All  equations  and  relationships  re¬ 
quired  to  predict  the  response  of  the  Turbine/Compressor  blade  to  a  spe¬ 
cified  loading  condition  were  compiled.  In  addition,  the  equations  re¬ 
quired  to  predict  the  amount  of  energy  dissipated  per  cycle  by  the  three 
forms  of  damping  were  developed.  The  experimental  model  was  used  to  gen¬ 
erate  data  for  several  variations  in  loading  conditions  for  the  purpose 
of  validating  the  model  for  use  as  a  too,  for  the  experimental  study  of 
slip  damping. 
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II.  Experimental  Development 


Model  Description 

Figure  1  is  a  picture  of  the  experimental  blade  model.  A  uniform 
beam  with  a  half  cylinder  fixed  to  one  end  was  chosen  as  a  simplified 
model  of  a  turbine  or  compressor  blade.  Due  to  the  limitations  of  this 
investigation,  and  the  phenomenon  being  studied,  this  simplification  was 
felt  to  be  acceptable.  The  half  cylinder  is  the  root  for  the  experimen¬ 
tal  blade  model.  The  shape  was  chosen  to  facilitate  both  the  manufac¬ 
turing  of  the  model  and  the  analytical  work.  The  dimensions  of  the  ex¬ 
perimental  model  are  shown  in  Figure  2.  Figure  3  is  a  picture  of  the 
experimental  disk  model.  The  assembled  blade-disk  model  (Figures  4  and 
5)  provides  a  simple  representation  of  a  turbine  or  compressor  blade 
mounted  in  an  engine  disk.  In  a  jet  engine  the  primary  force  acting  on 
the  blade  root  is  due  to  the  inertia  load  of  the  blade  on  the  blade  root 
which  is  produced  by  the  disk  and  blade  rotating  during  engine  operation. 
To  simulate  this  loading,  in  the  blade-disk  model,  a  spring  beam  with  a 
pointed  screw  was  used  to  apply  a  point  load  to  the  blade  root  at  the 
point  of  rotation  for  the  blade  root  (Figure  6) .  The  beam  was  instru¬ 
mented  on  both  sides  at  a  chosen  point  with  strain  gages  to  allow  for 
the  determir-acion  of  the  clamping  force.  The  blade-disk  model  was  manu¬ 
factured  using  Inconel-X.  This  material  was  chosen  because  it  is  repre¬ 
sentative  of  the  materials  used  in  a  jet  engine  and  to  allow  the  use  of 
the  blade-disk  model  at  elevated  temperatures  without  severe  changes  in 
the  material  properties.  The  properties  for  Inconel-X  are: 

6  2 

Modulus  of  elasticity  (E)  =  31  x  10  lbf/in 

Poisson's  ratio  (v)  =  .29 
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Figure  3.  Picture  of  Experimental  Disk  Model 
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Figure  4.  Picture  of  Assembled  Disk-Blade  Model 
(side  view) 


Figure  5.  Picture  of  Assembled  Disk-Blade  Model 
(end  view) 
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1.1 


Static  coefficient  of  friction  (m  )  = 

s 

Dynamic  coefficient  of  friction  =  -5 

Test  Set-up 

Figure  7  is  a  schematic  diagram  of  the  test  set-up  which  is  pictured 
in  Figure  8.  A  MB  model  N695  Automatic  Vibration  Exciter  Controller  was 
used  to  adjust  the  frequency  and  the  magnitude  of  the  signal  being  fed 
thru  the  MB  model  2120MB  Amplifier  to  the  MB  model  EA1500  Exciter.  The 
exciter  applied  a  sinusoidal  force  to  the  tip  of  the  blade  to  simulate 
the  dynamic  loading  of  a  compressor  or  turbine  blade  which  occurs  in  a 
jet  engine.  The  magnitude  of  the  applied  force  was  monitored  using  a 
Wilcoxon  model  Zll  impedance  head.  The  force  gage  signal  from  the  impe¬ 
dance  head  was  amplified  using  a  Kistler  model  568  charge  amplifier  and 
displayed  on  a  Tektronix  model  465M  dual  trace  oscilloscope.  Due  to 
problems  with  the  integral  accelerometer  within  the  impedance  head,  a 
separate  accelerometer,  mounted  on  the  blade  tip,  was  used  to  measure 
blade  tip  displacement.  The  signal  from  the  MB  model  306  accelerometer 
was  amplified  using  a  Kistler  model  504D  charge  amplifier  and  displayed 
on  the  same  scope  as  the  force  gage  signal.  The  signal  from  the  strain 
gages  was  monitored  using  a  Tektronix  type  549  oscilloscope  with  a  plug¬ 
in  dynamic  strain  preamplifier. 

Test  Procedure 

Due  to  the  nature  of  the  test  set-up,  with  the  blade  tip  constrained 
by  the  physical  attachment  of  the  exciter,  it  was  not  possible  to  study 
the  displacement  of  the  blade  tip  to  a  known  force  at  resonance.  Instead 
the  force  required  to  maintain  a  set  displacement  at  resonance  was  moni¬ 
tored.  The  procedure  for  collecting  the  data  was: 
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1)  The  force  gage  and  the  accelerometer  were  calibrated  using  the 
procedure  outlined  in  Appendix  B. 

2)  The  spring  beam  strain  gages  were  calibrated  using  the  procedure 
of  Appendix  C. 

3)  A  force  was  applied  to  the  blade  root  using  the  set  screw  and 
spring  beam. 

4)  A  displacement  for  a  data  run  was  selected. 

5)  A  frequency  was  selected  for  a  data  point  within  the  data  run. 

2 

6)  Using:  displacement  =  acceleration/ (2jt  frequency) 

the  acceleration  required  for  the  desired  displacement  at  the  selected 
frequency  was  determined. 

7)  Using  the  controller/oscillator,  the  needed  acceleration  was  set. 

8)  The  magnitude  of  the  force  being  exerted  to  maintain  the  displace¬ 
ment  was  recorded. 

9)  A  new  frequency  was  selected  and  steps  6,  7,  and  8  were  repeated. 
This  was  done  at  as  many  data  points  within  a  set  displacement  data  run 
as  desired. 

10)  Steps  4  thru  9  were  repeated  for  different  values  of  the  clamp¬ 
ing  pressure. 

Since  the  frequency  range  of  interest  was  the  range  about  a  point  of 
resonance,  the  frequency  range  for  the  data  runs  was  selected  about  the 
first  resonant  frequency  calculated  which  was  calculated  in  Appendix  A. 

Data  and  Discussion 

Figure  9  is  a  plot  of  the  force  versus  frequency  for  a  clamping  pres¬ 
sure  of  ten  pounds  and  three  variations  in  the  value  of  blade  tip  dis¬ 
placement.  The  first  item  to  note  is  the  presence  of  the  resonance  at 
320  cps .  This  corresponds  almost  exactly  to  the  value  of  resonance  which 


Frequency  (  CPS  ) 


Figure  9.  Plot  of  Force  Versus  Frequency  for  P  »  10  Ids  and 
Variations  in  Tip  Displacement 
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was  predicted  by  the  analytical  model  for  the  first  mode  of  the  blade. 

All  data  runs  were  made  over  the  range  of  240  to  420  cps,  but  a  sweep 
was  made  to  a  frequency  of  1200  cps  to  check  for  other  possible  reso¬ 
nances.  This  sweep  verified  the  existence  of  a  resonance  at  1000  cps, 
within  five  per  cent  of  the  value  of  resonance  which  was  predicted  by 
the  analytical  model  for  the  second  mode  of  the  blade. 

With  an  increase  in  the  blade  tip  displacement  there  is  an  increase 
in  the  magnitude  of  the  blade  tip  driving  force,  for  the  displacements 
of  .00005  in.  and  0001  in.  Figure  9  shows  an  amplification  factor  of 
1  to  6.5  at  320  cps.  This  indicates  no  change  in  the  damping  of  the  sys¬ 
tem.  This  is  difficult  to  understand  in  light  of  the  fact  that  all  the 

damping  mechanisms  are  functions  of  the  exciting  force.  The  slip  damp¬ 
ing,  however,  will  not  be  present  until  the  exciting  force  reaches  a 
threshold  value.  The  third  displacement  data  point  of  .0003  in.  plotted 
in  Figure  9  shows  an  amplification  factor  of  1  to  3.5  at  320  cps.  This 

indicates  a  significant  increase  in  the  amount  of  damping,  present  in 

the  system,  at  resonance.  It  may  be  assumed  that  the  threshold  value 
for  slip  damping  has  been  passed  and  that  slip  damping  is  making  a  large 
contribution  to  the  damping  of  the  system. 

Figure  10  is  a  plot  of  the  force  versus  frequency  for  a  clamping  pres¬ 
sure  of  twenty  pounds  and  three  variations  in  the  blade  tip  displacement. 
The  analytically  predicted  resonance  at  320  cps  is  present  in  this  set 
of  data.  The  amplification  factors,  for  the  resonance  peak  at  320  cps, 
for  the  lower  two  displacement  cases  are  almost  identical.  The  amplifi¬ 
cation  factor  for  the  largest  displacement  case  shows  a  significant  in¬ 
crease  in  the  amount  of  damping  present  at  resonsance.  Again  it  can  be 
assumed  that  the  critical  value  of  the  exciting  force  has  been  passed  and 
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that  slip  damping  is  contributing  to  the  dissipation  of  energy  of  the 
system. 

Figures  11  and  12  are  force  versus  frequency  plots  for  a  set  dis¬ 
placement  and  variations  in  the  clamping  force.  From  these  plots  it  can 
be  seen  that,  for  a  set  displacement,  or  exciting  force,  that  the  lower 
clamping  pressure  provides  a  better  reduction  of  the  amplification  fac¬ 
tor.  This  supports  the  conclusion  that  the  lower  the  clamping  pressure, 
the  more  slip  damping  is  present.  The  exception  to  this  occurs  at  320  cps 
on  Figure  11.  The  resonant  amplification  factors  are  almost  identical. 
This  can  be  explained  by  supposing  that  the  particular  value  of  the  ratio 
of  exciting  force  to  clamping  pressure  is  not  high  enough  for  slip  to 
occur.  The  prominant  resonance  at  370  cps  appears  to  correspond  to  the 
second  torsional  mode  of  the  blade. 
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III.  Analytical  Development 

Model  Description 

Figure  2  shows  the  dimensions  and  sign  convention  of  the  analytical 
blade  model  that  was  used  in  this  study.  The  blade  was  represented  as 
a  continuous  beam  of  constant  mass  per  unit  length  and  constant  cross 
sectional  area. 

Contact  Surface  State  of  Stress 

The  first  step  toward  the  calculation  of  the  energy  dissipated  by 
slip  is  the  calculation  of  the  stresses  acting  on  the  contact  surface. 

To  determine  the  traction  distribution  on  the  blade  contact  surface  re¬ 
quired  to  maintain  equilibrium,  the  loading  on  the  blade  root  was  broken 
into  three  cases  (see  Figure  13} .  Case  one  is  the  force  applied  to  the 
bottom  of  the  blade  root  by  the  setscrew.  This  force  was  intended  to 
simulate  the  inertia  force  of  a  blade  which  acts  on  a  turbine/compressor 
blade  root  due  to  the  rotating  motion  of  the  blade-disk  structure.  Cases 
two  and  three  are  the  transverse  force  and  moment,  respectively,  which 
are  applied  to  the  blade  root  due  to  the  dynamic  force  applied  to  the 
blade  tip.  The  first  two  loading  conditions  were  analyzed  using  the  solu¬ 
tions  for  the  stresses  on  a  wedge  loaded  at  it's  apex.  The  half  disk 
geometry  of  the  model  blade  root  can  be  matched  by  expanding  the  half 
angle  of  the  wedge  to  90  degrees.  For  the  third  loading  case  the  solution 
for  a  disk  loaded  by  a  couple  was  used.  By  starting  with  the  stress  func¬ 
tions,  for  the  wedge  and  disk  solutions,  and  evaluating  the  constants 
based  on  the  conditions  of  the  blade  root  model,  a  good  approximation 
of  the  stress  condition  of  the  contact  surface  can  be  obtained. 

Timoshenko  (23)  discusses  the  wedge  and  disk  problems  in  great  detail. 
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Figure  13.  Breakdown  of  Applied  Loads 
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1.  Case  One  Loading 


Using  as  the  stress  function 


4s  i  =  A^r6cos6 


(1) 


the  following  stresses  were  derived. 


o  =  -  A, sinfl 

r  r  1 


(2) 


UU  =  Tr0  =  ° 


To  evaluate  the  constant  A^,  the  integral  of  the  vertical  component  of 


stress  over  the  contact  surface  was  set  equal  to  the  applied  force  (P) 


/ 


r  0. 


o 


/ 


0 


/ 


a  r  sinQdO  +  la  r  sin0d0 
r  /  r 


TI-0. 


dz  =  -P 


(3) 


Substituting  equations  (2)  into  (3)  and  evaluating 
t 


2A 


J 


( *j0-l*sin20 ) 


+  (*30-f»sin20) 


0 


T1  -0  , 


dz  =  P 


(4) 


results  in  an  expression  for  A 


A.  = 


(5) 


1  t(201-sin2B1) 


Using  the  above  gives  the  stresses  on  the  contact  surface  due  to  case 
one  loading. 

2P  sinO 


o 


r  rt (sin2G^-20^) 


(6) 


°e  =  Tr0  =  0 


2.  Case  Two  Loading 
Using  as  the  stress  function, 
=  A„rOsin0 


'  2  '  2 

the  following  stresses  were  derived. 


(7) 


2 


0 


a  =  t 

0  re 


o  =  —  A„  cosG 
r  r  2 


(8) 


To  evaluate  the  constant  A^ •  the  integral  of  the  horizontal  component 
of  stress  over  the  contact  surface  was  set  equal  to  the  applied  force 
(F  sinojT)  . 


t  f  e 


a  r  co sudb  + 

r 


TI 

Tt-0  , 


r  cosede 


dz  =  -F  sinwT 


(9) 


results  in  an  expression  for  • 

. _ F  sinoiT  ,ln. 

2  t(2e^  +  sin2U^) 

Using  the  above  gives  the  stresses  on  the  contact  surface  due  to  case 
two  loading. 

2F  sinuiT  cosO 


°r  rt(26,  +  sin20,) 


(11) 


0 


=  Tr0  =  0 


3.  Case  Three  Loading 
Using  as  the  stress  function. 


>3  -  A38 


the  following  stresses  were  derived. 

o  =  a.  =  0 

r  0 


x  =  A3 
ro  ~ 

r 


(12) 


(13) 


Solving  for  by  setting  the  integral  of  the  shear  over  the  area  of 
the  contact  surface  equal  to  the  moment  caused  by  the  blade  tip  load 
(FL  sinwT  =  FM) . 
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Substituting  equations  (13)  into  (14)  and  evaluating, 


results  in  an  expression  for  . 

.  -FM 

A3  2t01 


(14) 


(15) 


(16) 


Putting  equation  (16)  into 
case  three  loading. 


°r  =  °e  =  0 

-FM 


re  2r261  t 


equations  (13)  gives  the  stresses  due  to 


(17) 


The  results  from  the  three  separate  loading  conditions  can  be 
combined,  using  the  principle  of  superposition,  to  provide  an  expres¬ 
sion  for  the  total  stress  condition  of  the  contact  surface. 


o  = 


2P  sinQ 


2F  sirnuT  cost) 
r  ”  rt(sin2Gi  -26^  +  rt (20^  +  sin2Gi) 

-FM 


r0  „  2  , 
2r  t6. 


(18) 


(19) 


°0  =  0 


(20) 


As  the  blade  loading  of  the  blade  root  increases,  the  shear  stress 
on  the  contact  surface  of  the  blade  root  increases  in  magnitude.  As  the 
value  of  the  shear  stress  at  a  point  on  the  contact  surface  readies  a 
value  equal  to  the  friction  force,  opposing  motion,  slip  will  occur. 

The  friction  force  is  generated  by  the  radial  stress  times  the  static 
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coefficient  of  friction.  Slip  will  continue  until  the  shear  force,  at 
the  point,  drops  to  a  value  which  is  less  than  the  friction  force  caused 
by  the  radial  stress  times  the  dynamic  coefficient  of  friction. 

A  non-dimensional  representation  of  the  stress  distribution  on  the 
contact  surface  is  pictured  in  Figure  14.  The  direction  of  the  exciting 
force  is  also  depicted  in  the  figure.  It  can  be  seen  from  the  figure 
that  the  normal  force  is  not  independent  of  the  exciting  force.  The 
static  and  dynamic  friction  forces  will  always  be  negative  while  the  shear 
force  will  change,  from  positive  to  negative,  depending  on  the  direction 
of  the  exciting  force. 

The  energy  dissipated  by  slip  is  the  integral,  over  the  slip  region, 
of  the  friction  force  times  the  relative  displacement  across  the  contact 
surface.  If  the  normal  force  was  independent  of  the  exciting  force,  then 
the  relative  displacement  across  the  contact  surface  could  be  related 
to  the  change  in  the  stress  component,  parallel  to  the  contact  surface, 
across  the  contact  surface.  Instead  the  relative  displacment  must  be 
calculated  by  first  solving  for  the  displacements  associated  with  the 
contact  surface  stresses. 


Contact  Surface  Displacements 

Associated  with  the  three  loading  cases  previously  discussed  are 
the  displacements  which  occur  at  points  on  the  contact  surface.  To 
solve  for  the  displacements  the  equations  for  the  relationships  of 
stress  to  strain,  and  strain  to  displacement  were  used. 


e 

r 


-vae) 


(21) 


e 


G 


U  1_  8V 

r  +  r  36 


E 


(a  -vo  ) 
6  r 


(22) 
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Direction  of  Loading  for 
Stress  Distribution  Shown 


Figure  1U .  Typical  Stress  Distribution  of  Contact 
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bra  r  ae  ar  r  G  rr0 
1.  Case  One  Loading 

Using  equations  (6)  and  equation  (21)  the  U  displacement  is 

u  =  -  123-JL.  +  f(u) 

tE(20  -  sin2G1)  '  1 


Using  equations  (24)  and  (22)  the  V  displacement  can  be  solved  for. 


_  2uP  cosO _  2P  cos 

tE (20  -  sin28  )  “  tE(20 


osO  log  r  /, 

;  -  w;,2V  -  y 


f(ii)dU  +  £(r) 


To  solve  for  f(0)  and  f(r),  equation  (23)  and  equation  (6)  are 
used  along  with  equations  (24)  and  (25) .  Using  these  equations  an  ex¬ 
pression  in  terms  of  f ( 0)  and  f(r)  is  developed. 


2  'P  cos  8 


df  (0)  1.  _  2P  cosO _  3 f  ( r )  _ 2 -'P  cos 8 

39  r  rth(20i  -  sin28  )  +  3r  +  rtE(20^  -  sin20^) 


■h 


For  this  equation  to  hold  the  sum  of  the  0  terms  must  be  equal  to 
zero  and  the  sum  of  the  r  terms  must  be  equal  to  zero.  Timoshenko  (23) 
discusses  the  solution  to  equation  (26)  for  f(B)  and  f(r). 

f ( r)  =  Qjr  (27) 


f(B)  =  f (particular )  +  f (homogeneous)  (28) 
The  particular  solution  is  dependent  on  the  other  0  terms  present  in 
the  expression  for  f(o)  and  f(r),  equation  (26). 

f<0)  -  '  +  fi2fllnU  +  e3cos0  (29) 
Using  equations  (27)  and  (28) ,  the  displacements  are 


2P  sin'J  log  r  P  ( v  -  1)8  cosO  „  .  n  .  _  ,, 

tJ  tE(201  -  s i n 2 0 ^ )  “  tE(201  -  sin201)  +  C2Sin  Q3COfa  (30) 
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and 


2yP  cos6 _ 

tE(20  -  sin201) 


2P  cosu  log  r  P(y  -  1)  cosQ 

tE(201  -  sin201)  tE(201  -  sin20  > 


P(v  -  1)  sine 
tE(201  -  ain201) 


+  Q^cosQ  +  Q^sint)  +  Q^r 


(31) 


Selecting  for  a  boundary  condition  that  the  displacements  (U  and  VJ  are 
equal  to  zero  at  the  point  on  the  contact  surface  where  tt  =  0  allows  the 
evaluation  of  the  constants  in  equations  (30)  and  (31). 


°1  =  ° 

_  P  ( y  +  1) _  2P  log  r _ 

"  tE(26.  -  Sin261>  +  tE(26  -  sin2u.) 


(32) 
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Substituting  equations  (32)  into  equations  (30)  and  (31) 

P(v  -  1)0  cos(J  P(v  1)  sinO 

tE(201  -  sin20  )  tE(20  -  sin20  ) 


and 

P(v  -1) 0  sinO 
tE(20  -  sin201) 


(34) 


2.  Case  Two  Loading 

Using  equations  (11)  and  equation  (21)  the  U  displacement  is 
2F  sinooT  cosQ 


U  = 


tE(20  +  sin261) 


log  r  +  f ( 0 ) 


(35) 


Using  equations  (35)  and  (22)  the  V  displacement  is 

2F  sinwT  sin8  _  2vF  sintoT  sin8 

V  tE(20  +  sin201)  109  *  ~  tE(26  +  sin201) 


d0  +  f (r) 


(36) 


To  solve  for  f(0)  and  f(r)  the  same  method  is  employed  as  was  used 
for  the  case  one  loading  displacement  analysis.  The  results  are 

f ( r)  =  Q  r  (37) 
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and 


f  (0) 


(v  -  1)  F  sinmT  cosu 
tE(20^  +  sin2(i^) 


+  Q^sinG  +  Q-jcos0 


(38) 


These  results,  when  used  with  equations  (35)  and  (36)  provide  the 

following  displacement  expressions, 

2F  sinuiT  log  r  cost)  (v  -  1)  F  sinujT  cosu 
U  tE ( 20  t  sin26  )  +  tE{2ei  +  sin2ei> 


(39) 


+  Q„sin0  +  Q,cost) 
2  3 


and 


v  _  2F  sinwT  log  r  sinO  _  (v  -  1)  F  sinmT  cost) 
”tE(2&1  +  sin2«  )  ”  tE(2«1  +  sin2ei> 


(v  -  1)  F  sinmT  sine 


(40) 


(v)  "  r  aiimii  saho  .  .  „  _  ,s  _ 

tE(2o1  +  sin201)  °2COS0  "  ^3slne  +  Cir 
To  solve  for  the  constants  in  equations  (39)  and  (40)  the  displacement 
at  the  point  0  =  0  is  set  equal  to  zero.  Using  this  boundary  condition 
the  constants  are 


21  =  ° 

(v  -  1)  F  simoT 
y2  ~  tE(2i)  +  sin201) 

_  2F  sinmT  log  r 
2 3  tE(201  +  sin29^) 

Introducing  equations  (41)  into  equations  (39)  and  (40) 

( v  +  1)  F  sing)T  cos0  _  (y  -  1)  F  sinuT  sinO 
U  ~  tE(201  +  sin2G^)  “  tE(20i  +  sin201> 


and 


2  ( v  ~  1)  F  sinmT  cos6  (v  -  3)  F  sinmT  sinO 
tE(201  +  sin201)  ”  tEUl^  +  sin2«1) 


(41) 


(42) 


(43) 


3.  Case  Three  Loading 

Using  equations  (17)  and  equations  (21)  and  (22)  it  is  possible 
to  establish  the  relationships 
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U  =  f  (  0 ) 


V  =  -  I  f  (0)  do  (45 
Using  equation  (23)  it  is  possible  to  set  up  an  expression  in  terms 
of  f  (U) 


FM 

2G  r  t  0 , 


Of  (0) 


+  I  f(0)  d0 


which  can  be  solved  for  f(0) 


f(6)  =  Q2sin0  +  Q3cos0  -  e- 


where 


2(1  +  v) 

which  gives  the  following  displacements 


U  =  Q2sine  +  Q3cose  -  — 


V  =  Q2cos0  -  Q3sin0  *  '2q' 7"tT 


The  displacement  at  0  =  0  is  set  equal  to  zero  in  order  to  solve  for 
the  unknown  constants.  Solving  for  the  constants  Q and  Q3  and 
putting  them  into  equations  (49)  and  (50)  gives  the  final 
displacement  equations 

U  =  -^—7—  (cosO  -  1)  (51) 

E  r  t  03 


V  =  'e  rlT^  (U  '  Slnt)) 

The  total  displacements  can  be  arrived  at,  in  the  same  manner  as 
the  total  stress  components  were  achieved,  by  the  principle  of 


superposition. 


U 


and 


P(v  ~  1)0  cosq  Pjy  +  I )  sinO  (y  +  1)  F  sin^TE)  cos6 

tE (20  -  sin2o  )  tE<20,  ~  sin2b.)  tE(20,  +  sin20,) 

1  1  1  1  11 


(v  -  1)  F  sinujT  si_nlj  ( 1  +  v )  FM  ( cos 9  -  1) 
tE(201  +  sin2l>  )  E  r  t  6 x 


(53) 


P(v  -  1)  u  sin0  2  (v  -  1)  F  sinoT  cos0  (v  -  1)  F  sin^'TQ  sin^ 

tE  (20  x  -  sin20  “  tE(20  +  s^1^)  ”  tE(2Q;L  +  sin201) 

.  (1  +  v)  FM  (6  -  sinO)  lr.s 

+  - 7, -  (54) 

E  r  t  0 

Having  expressions  for  the  0 ,  or  V,  and  the  radial,  or  U,  displace¬ 
ments  it  is  possible  to  solve  for  the  relative  displacement  between  two 
originally  adjacent  points  after  slip  has  occured.  Figure  15  depicts 
the  quantity  that  is  desired. 

The  component  of  the  displacement  of  interest  is  the  0  direction  or 
V  component.  On  one-half  of  the  slip  region,  region  1,  the  displacemnt 
will  be  described  by  equation  (54).  Denoting  equation  (54)  in  that  half 
as 


V1  =  f(°r> 


+  g(Tr0> 


(55) 


where 


and 


f(  )  is  the  V  displacement  associated  with  the  radial  stress 


g(  )  is  the  V  displacement  associated  with  the  shear  stress 
In  the  other  half  of  the  slip  region,  region  2,  the  displacement  (V  ) 
will  be  a  combination  of  the  displacement  associated  with  the  radial 
stress  (f(  ))  and  the  displacement  associated  with  the  friction  force 
acting  on  that  portion  of  the  slip  region  (h(  )).  This  can  be  written 


as 


Figure 
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15.  Relative  Displacement  Across  Contact 
Region 
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V2  =  f(ar)  +  h(uD°r) 


(56) 


The  relative  displacement  is  the  difference  between  equations  (55)  and 
(56) 


V1  -  V2  '  9  '’rt1  -  h<Vr> 


which  can  be  written  as 


(57) 


V  -  V 
1  2 


(1  +  v)  FM  (6  -  sin1-))  lD  P(V  “  1)0  Sln6 


E  r  t  6. 


tE(201  -  sin261) 


(58) 


2(v  -  1)  E'  sinuiT  cos6  (y  -  1)  F  singiTO  sin8 
tE(2eL  +  sin201)  +  tE(2t)  +  sin281) 

Since  the  energy  dissipated  is  related  to  the  integral,  over  the 
slip  region,  of  the  friction  force  times  the  relative  displacement  it  is 
necessary  to  determine  the  limits  of  integration. 


Points  of  Slip  Initiation  and  Zero  Normal  Force 

There  will  be  two  such  points  (Q  and  R) ,  one  on  each  side  of  the 
blade  root,  where  slip  will  begin  (Figure  14) .  The  equation  for  the  two 
points  where  the  friction  force  and  the  shear  stress  are  equal  is 


-  ABS  1  =  no 

rQ  S  r 


(59) 


Solving  for  points  Q  and  R,  the  points  of  equality,  on  the  blade 
contact  surface  by  writing 


n_o  =  LsinB  +  Gcoso 
S  r 


(60) 


where 


2PU 


L  = 


rt(sin2h^  -  20^) 


(61) 


G  = 


2FUg  sinmT 
rt(sin20^  +  20^) 


and  writing 
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"  ABS  TrO  =  H  (62) 

setting  the  two  functions  equal  to  each  other 

LSinO  +  GcosO  =  H  • 
i 

rearranging  the  equaticfp  and  squaring  both  aides 

f 

L2 (1  -  cos20)  =  G2cos2U  +  H2  -  2HGcosu 
collecting  terms  results  in  a  quadratic  equation 


(L2  +  G2)  cos26  -  2HGcos6  +  (H2  -  L2)  =  0  (63) 

This  equation  can  be  solved  for  the  angles  corresponging  to  points 
Q  and  R.  1 

Due  to  the  nature  of  the  loading  on  the  blade  root  there  will  be  a 
small  portion  of  the  contact  surface  adjacent  to  the  y-axis  which  will 
be  free  of  any  normal  force.  This  small  portion  of  the  contact  surface 
will  have  an  outward  normal  opposite  in  direction  to  the  direction  of 


the  tangential  force  discussed  as  the  case  two  loading  condition.  To 
solve  for  the  point  (S)  where  the  normal  force  equals  zero, 

_ 2P  sinO  2F  sinurr  cos© 


rt(sin20^  -  26^) 


rt(20^  +  sin28^) 


(64) 


A 


Solving  the  equality  for  the  angle  of  the  point  (S) . 
0  =  tan'1  F<sin26i  "  2V 


P(20  +  sin281) 


(65) 


Having  an  expression  for  the  normal  force  on  the  contact  surface 
(equation  18)  the  friction  force  can  be  written  as 


M 


0 

D  r 


2Pud  sin0  2FUD  sinufT  cos6 

tr(sin20^  -  26^)  +  tr(sin20^  -  20^) 


(66) 


The  relative  displacement  across  the  contact  surface  due  to  slip  is 
known  (equation  58) ,  and  the  points  of  slip  initiation  are  known 


3? 


(equation  65) .  With  these  expressions  the  energy  dissipated  due  to  slip 
can  be  calculated. 


Energy  Dissipated  by  Slip 

The  energy  dissipated  by  the  slip  region  is  the  product  of  the  fric¬ 
tion  force  acting  over  the  region  times  the  relative  displacement  across 
the  interface.  This  quantity,  integrated  over  the  slip  region,  will  give 
the  amount  of  energy  dissipated  at  a  particular  instant.  Integrating  this 
result  over  one  cycle  will  result  in  the  energy  dissipated  per  cycle. 
During  one  cycle  the  exciting  force  will  change  from  zero  to  a  maximum 

value  four  times.  Thus  the  energy  dissipated  per  cycle  (D  ,  lb-in)  is 

o 


D  =  4 
o 


/ 


Friction  force  x  Relative  displacement 


(67) 


slip  region 

Four  times  the  product  of  equations  (66)  and  (58)  integrated  over  the 
slip  region  gives  the  energy  dissipated  per  cycle. 

Q  71 


D  =  4 
o 


/ 


....  y 


<  VV  <VX  -  V2)  d6 


(68) 


It  should  be  noted  that  the  region  of  zero  normal  force  is  included 
in  the  above  expression.  The  region  is  on  the  order  of  two  degrees  and 
its  inclusion  in  the  expression  has  a  negligible  effect. 

In  order  to  isolate  the  effect  of  the  slip  damping  it  is  necessary 
to  develop  expressions  for  the  air  damping  and  the  material  damping. 


Energy  Dissipated  Ry  Air  Damping 

The  amplitudes  associated  with  the  blade  displacements  are  small 
and  thus  the  force  acting  on  the  blade  is  proportional  to  velocity 
(reference  2) 
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To  develop  an  expression  for  the  hysteresis  damping  of  the  blade, 
a  relationship  between  stress  and  energy  dissipated  of  the  form 


D  =  J  CJ 


was  assumed.  J  and  N  are  empirically  derived  constants  and  is  the 

bending  stress  due  to  the  tip  load  on  the  blade.  D  has  units  of 
3 

(lb-in/in  )  and  must  be  integrated  over  the  total  valume  of  the  blade 
to  obtain  the  energy  dissipated 


t  b/2 


//  /<¥? 


dxdydE 


0  -b/2  0 
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Evaluating  this  integral  where  M  =  F(L  -  X)  leads  to 


D 

o 


2Jt 

(N  +  1) 


y 

N 

bL’ 

~2 

i 

2 

For  Inconel-X  and  less  than  33  KSI 


(75) 


N  =  2.224 

-13  (76 

J  =  4.586  x  10 

In  order  to  be  able  to  use  the  experimental  model  to  verify  the 
analytical  model,  it  is  necessary  to  develop  an  expression  which  re¬ 
lates  such  quantities  as  blade  tip  displacement,  exciting  force, 
clamping  pressure,  material  properties,  etc.  This  can  be  done  by 
developing  a  solution  to  the  equation  of  motion  for  the  blade.  The 
first  step  toward  this  is  the  evaluation  of  the  restoring  moment  due 
to  slip  which  acts  on  the  blade  root. 


Restoring  Moment  Due  to  Slip 

The  restoring  moment  due  to  slip  (M)  can  be  calculated  from 


M  =  r  f  I  /pa  d  0  +  /  p  o  d0 

J  [J  °r  J  °r  J 


dz 


(77) 


where  t  is  the  depth  of  the  blade  root  in  the  z  direction.  Placing 
equation  (18)  into  equation  (69)  and  evaluating, 

Q 


M  =  ru 


(sin20 


-2P  r 
i  -  2v  r 

2F  sininT _  I 

20^  +  sin20^)  I 


+  cosQ 


'R 


+  sin0 


;] 


This  leads  to  an  expression  for  the  restoring  moment, 


2P(cosR  -  cosQ 
M  ~  r,D  (sin201  -  201 


-  2)  2F  siniiJT 
)  +  (20 


(sinQ  -  sin  R) 
+  sin20^) 


(78) 


which  can  be  rewritten  as 


M  =  m  sinuiT 


where 

m 


2P(cosR  -  cosQ  -  2) 
(sin20^  -  20^)sin(ijT 


2F(sinQ  -  sinR) 
(201  +  sin261) 


(79) 


It  should  be  noted  that  the  region  of  zero  normal  force  was  in¬ 


cluded  in  the  determination  of  the  restoring  moment  which  was  due  to 


friction  developed  on  the  contact  surface.  The  region  is  on  the  order 
of  two  degrees  and  its  inclusion  in  the  integration  has  a  negligible 


effect  on  the  result. 


Equation  of  Motion  and  Boundary  Conditions 

To  derive  the  equation  of  motion  and  boundary  conditions  for  the 
blade,  Hamilton’s  principle  was  used. 

r 

J  (<5K  -  6U  +  6WK)dT  =  0  (80) 

0 


For  the  blade  the  variation  of  the  kinetic  energy  is 


6K  = 


f 


aw 

8T 


1 

3T 


~  <5W  dx  +  I 


a _ w 

o  3x3T 


x=0 


3x3T 


6W 


x=0 


(81) 


where 


L  is  the  length  of  the  blade  in  inches 

y  is  the  mass  per  unit  length  of  the  blade  (lbm/in) 

I  is  the  mass  moment  of  inertia  of  the  blade  root  about 
o 

2 

point  0  (lbm-in  ) 

W  =  W(x,T)  is  the  displacement  of  the  blade  in  the  Y  direction 
The  variation  of  the  potential  energy  of  the  blade  is 


fiU  = 


a2w  a2 

„  2  2 
dx  3x 


6W  dx 


(82) 
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where 


2 

E  is  the  modulus  of  elasticity  of  the  blade  (lbf/in  = 
386  lbm/sec-in) 

4 

I  is  the  area  moment  of  inertia  of  the  blade  (in  ) 
and  the  variation  in  the  work  is 


6WK  =  F  sinwT  tW 


-  m  sinwT  t —  6W 
x=L  3X  x=0 


(83) 


Placing  equations  (81),  (82),  and  (83)  into  equation  (80),  performing 

the  necessary  integrations,  and  eliminating  the  zero  terms  leads  to 
the  equation  of  motion. 


4  2 

aw  J  aw 

EI  ~  +  y  — -  =  0 

3x  3T 


the  boundary  conditions  at  x=0, 
W  (0  ,T)  =  0 

and 


EX  —  W  (0  ,  T) 

ax 


—  W  ( 0  , T )  II  +  m  sintoT  =  0 


and  the  boundary  conditions  at  x  =  L 

,2 

EI  — — —  W  (L,T)  =  0 
3x 


and 


EI  W(L,T)  +  F  sinoiT  =  0 
3x3 


(84) 


(85) 


(86) 


(87) 


(88) 


Solution  to  Equation  of  Motion 

To  solve  the  equation  of  motion  for  the  displacement  it  is  neces¬ 
sary  to  make  the  boundary  conditions  homogeneous.  Denoting 


3W 

3x 


Vi 

x 


W 

XX 


etc. 
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and  writing  the  equation  of  motion 


4 

W  -  a  V* 
xxxx  tt 


0 


where 


Making  a  change  of  variables 

W(x,T)  =  Y(x,T)  +  (B^sinaH'X  +  B^cosa'J'x  +  B^sinhavx 

2 

+  B  ,  coshail'x)  sin'i'  T 
4 

with 

i|/2  =  w 

From  equation  (91) ,  and  the  boundary  condition  (85)  for 
W(x,T)  =  Y(x,T)  =  0 


Placing  the  result  in  equation  (91)  and  substituting  this  into 
(87)  and  requiring  that 


W  (L,T)  =  Y  (L,T)  =  0 
XX  XX 

requires  that 

B^sinha'/'L  =  B^sina^L  +  B 2  (cosa^L  +  cosha'I'L) 
Using  equation  (88)  and  forcing 


W  (L,T)  =  Y  (L,T)  =  0 
xxx  xxx 


requires  that 


B^osa^L  -  B  (sina^L  -  sinha'i'L)  -  B  cosha'I'L  =  3.3 

1  z  J  fc-la  v 


Using  equation  (86)  and  forcing 

El  W  (0,T)  +  I  Wxtt (0 , T)  =  El  Y  (0,T)  +1  W  (0,T)  =  0 
xx  o  xx  o  xtt 

requires  that 

m  -  I  ail5<0,  +B  i 
„  o  1  3 


(89) 


(90) 


(91) 


equation 


(92) 


(93) 


(94) 
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Equations  (92),  (93),  and  (94)  provide  three  equations  in  three 
unknowns  which  can  be  solved  to  provide  the  following  results. 


B1  = 


m 

F 

^  (  sinai^L 

- 

T  5 

1  a 

3  3 

Ela  ^  (cosai^L  -  cothayL  sina^L) 

sinha^L 

o 

. 

c 


,  s  i  nha'IL 

1  + : - t —  + 

sinaVL 


[ 


1  + 


sinal’L 


cosa'lL  +  coshaij'L 


(sinaiR,  -  sinha’^I,)  -  (cosayL  +  ooshaVL)  cothatlL 


cosa'IL  -  cothaij'L  sinailL 
F 


B2  = 


EIa^4’^  (cosaj'L  -  cothaij’L  sina^D 
F 


(95) 


I  (|^a 
o 


Lr 

x 

r 


3  3 


E 


B3  = 


L 

«■ 

1  (  sinhailL 

sinaiJ'L  ( 

1  sina^L 

cosa^L  +  coshatJ'L 

+  sinailL  I 
sinha^L I 

J 

i 


Ela  <l>  (cosa<i<L  -  cotha^L 

_ cosajh  -  cotha'I'h  sina <PL 

sinatH,  -  sinhailL)  -  (cosa'vL  +  cosha^L)  cotha'lLj 

(96) 


(cosa^L  ~  cothailL  si  na<lL)  sinhailb 
(cosal<L  +  cosha'IL) 


sinha^L 


B. 


F  sinaij'I, 


(97) 


3,  3 


Ela  if*  (cosa^L  -  cothaiJ  L  sina^L)  sinha^L 
With  the  change  in  variables,  and  the  constants  (95) ,  (96) ,  and  (97) 


the  homogeneous  equation  of  motion  is 

Y  +  Y  =0 
xxxx  El  tt 

and  the  homogeneous  boundary  conditions  are 

Y  (0 ,  T)  =  0 


Y  ( L,T)  =  0 
xx 


Y  (L,T)  =  0 

XXX 


il  Y  ( 0 , T )  +  I  Y  (0,T)  =  0 

XX  o  xtt 


(98) 

(99) 
(100) 

(101) 

(102) 
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Assuming  a  solution  of  equation  (98)  of  the  form 
Y(x,T)  =  X(x)Z(T) 

and  using  the  solution  method  of  separation  of  variable  results  in 

Z(T)  =  JcoswT  +  KsinuiT  (103) 

and 


X(x)  =  D^sinfix  +  D^cosjsx  +  D^sinhyx  +  D^coshhx 
where 


(104) 


(105) 


Using  the  first  boundary  condition  (99) 


°2  =  '°4 

Using  the  second  boundary  condition  (100) 

-  D^singL  +  (cosgL  +  cosh(3L)  +  D^sinhgL  =  0 

Using  the  third  boundary  condition  (101) 

-  D_  cosfcSL  +  DfsinhtfL  -  sin0L)  +  D^coshBL  =  0 

14  3 

Using  the  fourth  boundary  condition  (102),  rewritten  as 

I  (i2 
o 


(106) 


(107) 


(108) 


X  (0) - —  X  (0)  =  0 

xx  El  x 


gives 


I  u 

2d4  -  -fie  <Di  *  d3> 


(109) 


For  the  non-trivial  solution  where  D  ,  D^,  D^,  and  are  not  zero, 
the  determinant  of  the  coefficient  matrix  must  be  zero.  The  governing 
equation  is 

I  u,2 

— -  (cosht>L  cosBL  t  1)  +  cosBb  sinhBL  -  sintfL  coshBb  =  0  (110) 

El  P 

This  equation  can  be  solved,  to  provide  the  values  of  tLandio,  using 

S 


( BL) 4EI 
YL4 


(111) 


40 


I 


See  Appendix  A  for  the  solutions  to  (110)  and  (111).  Using  equations 

(107) ,  (108) ,  and  (109) ,  the  constants  D.,  and  D„  can  be  put  in  terms 

3  4 

of  to  provide 


D3  =  D, 


3q | (costjL  +  coshBL)  +  sinBL 


2L3y 


+  -2±~±- 


sinhfcSL 


sinBL 


sinhBL 


3  (cosBL  +  coshBL) 

I  v 
O 


(112) 


and 


°4  =  °1 


*F^] 


(cosBL  +  coshgL) 
sinhBL 


1  (BL)  .  OT 

o  sinBL 

t  - 1 — : — ; -  + 


sinBL 


3  sinhgL  (cosBL  +  coshBL) 
Z.L  y 


(113) 


Setting  =  1  and  substituting  (112),  (113),  and  (106)  into  equation 
(104)  provides. 


XN(x)  =  sinBNx  + 


1  ( Bl)  sin ( Bl) 
+  o  N  _ N 


sin  (3  L) 

N 


3  sinh  ( 8  L)  (cos  (6  L)  +  cosh  (3D) 
2L  Y  N  N  N 


x  (coshB  ,x  -  cosB  x) 
N  N 


I  (  3D  (cos  (3D  +  cosh  (3  L) ) 
o  N  N  N 


2L3  y 


2L3  Y 


sinh ( B^L) 


sin  ( 3D 
N 


T  ti3  (cos  ( 8  L)  +  cosh  (3D) 
I  ( BL)  N  N 

O  N 


sin ( 3 L) 

sinh  ( 3,,L) 
N 


(sinh 3  x) 
N 


(114) 


Using  equations  (114)  and  (103),  the  solution  to  equation  (98)  can  be 


2 

1 

Io<V') 

(cos  ( 3..L)  +  cosh  (3  L))  1 

N  N  - 

4 

YL3 

sinh  (3D 

N 

. 

41 


1 


written 


Y(x,T) 


(J  cosiLT  +  K  sinfiT) 
N  N 


(115) 


Substituting  the  above  into  equation  (91)  provides  an  expression  for 
the  displacement. 


W  ( x ,  T )  =  X  (x)  (J  cosuiT  +  K  sinu  T) 

N  N  N  N  N 

N=1 

+  ^B^sinatyx  +  (cosaiJ.'X  -  coshaijix) 

+  B^sinhaiJjx  J  sinfi^T)  (116) 

To  evaluate  the  remaining  constants  J  and  K  it  is  necessary  to  use 

N  N 

the  initial  conditions 

W(x,0)  =  0  (117) 

and 

wfc(x,0)  =  0  (118) 

The  first  initial  condition  (117)  requires  =  0.  Using  this  result 
and  the  second  initial  condition  (118)  requires  that 

2 

K  <3  X  (x)  =  -  ip  (B  Sinai Jix  +  B  (cosaj-x  -  cosha^x) 

N  N  N  1  Z 

+  B^sinhai|ix)  (119) 

The  functions  <h  X  are  not  orthogonal  functions.  Therefore  a 
N  N 

collocation  method  must  be  used  to  solve  for  the  constants  (1C,)  . 

N 

The  desired  number  of  modes  (N)  must  be  chosen,  and  the  X  (x)  as 

N 

well  as  the  right  hand  function  must  Vie  evaluated  at  a  number  of  points 
along  the  blade.  The  X^(x)  can  be  evaluated  using  equation  (110)  and 
the  right  hand  function  is  evaluated  using  equations  (95),  (96),  and 

(97).  A  matrix  expression  is  developed 
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constants 


(120) 


fj(X) 

V*> 

which  can  be  solved  for  the  K  's. 

N 

The  solution  to  the  equation  of  mot  ion  provides  an  expression  which 
relates  the  various  parameters  of  the  experimental  and  analytical  models. 
This  expression  can  be  used  as  a  moans  of  verifying  the  analytical  model 
with  the  experimental  model.  It  can  also  he  used  study  the  effects 
of  varying  the  different  parameters.  If  there  is  or re lation  between 
the  analytical  model  and  the  experimental  model,  then  the  expressions 
for  the  energy  dissipated  can  be  assumed  to  be  correct . 
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IV.  Recommendations 


The  recommendations  are  put  forth  in  hopes  that  some  useful  infor¬ 
mation,  regarding  slip  damping  of  Turbine/Compressor  blades,  can  be 
generated,  using  the  analytical  and  experimental  models  provided  by 
this  study. 

1)  Using  the  analytical  model  the  following  quantities  for  various 
combinations  of  exciting  force  and  clamping  pressure  should  be  generated 

a)  Value  of  exciting  force  for  complete  slip. 

b)  The  exciting  force  necessary  to  obtain  a  set  blade  tip  dis¬ 
placement  . 

c)  A  comparison  of  the  energy  being  dissipated  by  the  three 
forms  of  damping. 

2)  Using  the  experimental  model 

a)  A  wider  variation  of  clamping  pressures  and  exciting  force. 

b)  Differing  conditions  for  the  contact  surface,  such  as  a 
sandwiched  viscoelastic  layer. 

c)  Investigations  of  slip  in  an  environment  with  an  elevated 
temperature . 

3)  The  analytical  model  should  be  modified  to  include  a  sandwiched 
viscoelastic  layer  at  the  contact  surface. 
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Appendix  A 


Solution  To  Characteristic  Equation 


The  characteristic  equation 


I  *2 
o 

El 


(coshBL  cosBL  +  1)  +  cosBL  sinh&L 

-  sinBL  coshBL  =  0 


(A-l) 


can  be  solved  using  a  combination  of  the  bisection  and  secant  methods 
to  find  the  points  (B^L)  where  equation  (A-l)  is  satisfied.  Using 
the  relation 


„2  B  El 

w  =  - 

Y 


and  the  values 

L3  =  1157.625  in3 

I  =  .01828  lbm-in2 
o 


(A-2) 


Y  =  .18625 


lbm 

in 


The  first  five  values  of  BL  which  satisfy  equation  (A-l)  are  listed 
in  column  two  of  Table  One.  In  column  three  of  the  same  table  are 
listed  the  classic  values  of  BL  for  a  clamped-pinned  beam.  Using 
equation  (A-2)  the  frequencies  corresponding  to  the  different  values 
of  BL  are  entered  in  column  four  of  Table  One. 
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Appendix  B 


Accelerometer  and  Impedance  Head  Calibration 

Figure  16  shows  the  set-up  used  to  calibrate  the  accelerometer 
and  force  gage.  The  large  MB  model  CIO  exciter  and  the  MB  model  3402MB 
amplifier  were  used.  The  accelerometer  and  the  force  gage  were  attached 
to  the  head  of  the  exciter.  A  known  mass  was  placed  over  the  impedance 
head.  The  signals  from  the  accelerometer  and  the  force  gage  were  fed 
thru  the  same  data  chain  as  was  used  for  the  data  taking.  A  known 
acceleration  was  exerted  by  the  exciter.  The  output  of  the  accelero¬ 
meter  and  the  force  gage  were  recorded.  Knowing  the  mass  over  the  force 
gage  allowed  the  calculation  of  the  force  corresponding  to  a  certain 
signal  level. 
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Figure  16.  Schematic  of  Calibration  Set-up  for  the 

Accelerometer  and  Impedance  Head  Calibration 
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Appendix  C 


To  calibrate  the  spring  beam  it  was  necessary  to  apply  a  known 
force  and  record  the  corresponding  signal  from  the  strain  gages. 

Figure  17  shows  the  force  table  and  guide  which,  when  assembled  on  the 
disk  model  (Figure  18) ,  allowed  a  known  force  to  be  applied  to  the 
spring  beam.  Increasing  amounts  of  weight  were  set  on  the  force  table 
and  the  signal  from  the  strain  gages  was  recorded. 


figure  17.  Picture  of  Force  Table  and  Guide 


Figure  18.  Picture  of  Force  Table  in  Place 
on  the  Disk  Model 
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An  experimental,  and  an  analytical,  model  of  Turbine/Compressor  blade 
with  a  circular  root  geometry  is  developed  for  the  purpose  of  studyinq 
the  effects  of  slip  damping  on  the  response  of  the  blade.  The  expression 
for  the  energy  dissipated  by  slip  damping  is  presented.  In  addition  the 
expressions  for  the  energy  dissipated  by  air  and  hysteresis  damping  are 
presented.  The  equation  of  motion  for  the  blade  is  developed,  and  solved. 
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to  provide  an  expression  relating  the  various  parameters  of  the 
experimental  and  analytical  models. 


